Synthesizing environmental, epidemiological and vector and parasite genetic data to assist decision making for disease elimination.

We present a framework for identifying when conditions are favourable for transmission of vector-borne diseases between communities by incorporating predicted disease prevalence mapping with landscape analysis of sociological, environmental and host/parasite genetic data. We explored the relationship between environmental features and gene flow of a filarial parasite of humans, Onchocerca volvulus, and its vector, blackflies in the genus Simulium. We generated a baseline microfilarial prevalence map from point estimates from 47 locations in the ecological transition separating the savannah and forest in Ghana, where transmission of O. volvulus persists despite onchocerciasis control efforts. We generated movement suitability maps based on environmental correlates with mitochondrial population structure of 164 parasites from 15 communities and 93 vectors from only four sampling sites, and compared these to the baseline prevalence map. Parasite genetic distance between sampling locations was significantly associated with elevation (r = .793, p = .005) and soil moisture (r = .507, p = .002), while vector genetic distance was associated with soil moisture (r = .788, p = .0417) and precipitation (r = .835, p = .0417). The correlation between baseline prevalence and parasite resistance surface maps was stronger than that between prevalence and vector resistance surface maps. The centre of the study area had high prevalence and suitability for parasite and vector gene flow, potentially contributing to persistent transmission and suggesting the importance of re-evaluating transmission zone boundaries. With suitably dense sampling, this framework can help delineate transmission zones for onchocerciasis and would be translatable to other vector-borne diseases.

Cytotaxonomic characterization and estimation of migration patterns of onchocerciasis vectors (Simulium damnosum sensu lato) in northwestern Ethiopia based on RADSeq data.

BACKGROUND: While much progress has been made in the control and elimination of onchocerciasis across Africa, the extent to which vector migration might confound progress towards elimination or result in re-establishment of endemism in areas where transmission has been eliminated remains unclear. In Northern Ethiopia, Metema and Metekel-two foci located near the Sudan border-exhibit continuing transmission. While progress towards elimination has been faster in Metema, there remains a problematic hotspot of transmission. Whether migration from Metekel contributes to this is currently unknown. METHODOLOGY/PRINCIPLE FINDINGS: To assess the role of vector migration from Metekel into Metema, we present a population genomics study of 151 adult female vectors using 47,638 RADseq markers and mtDNA CoI sequencing. From additional cytotaxonomy data we identified a new cytoform in Metema, closely related to S. damnosum s.str, here called the Gondar form. RADseq data strongly indicate the existence of two distinctly differentiated clusters within S. damnosum s.l.: one genotypic cluster found only in Metema, and the second found predominantly in Metekel. Because blackflies from both clusters were found in sympatry (in all four collection sites in Metema), but hybrid genotypes were not detected, there may be reproductive barriers preventing interbreeding. The dominant genotype in Metema was not found in Metekel while the dominant genotype in Metekel was found in Metema, indicating that (at the time of sampling) migration is primarily unidirectional, with flies moving from Metekel to Metema. There was strong differentiation between clusters but little genetic differentiation within clusters, suggesting migration and gene flow of flies within the same genetic cluster are sufficient to prevent genetic divergence between sites. CONCLUSIONS/SIGNIFICANCE: Our results confirm that Metekel and Metema represent different transmission foci, but also indicate a northward movement of vectors between foci that may have epidemiological importance, although its significance requires further study.

Genome-based tools for onchocerciasis elimination: utility of the mitochondrial genome for delineating Onchocerca volvulus transmission zones.

National programs in Africa have expanded their objectives from control of onchocerciasis (river blindness) as a public health problem to elimination of parasite transmission, motivated by the reduction of Onchocerca volvulus infection prevalence in many African meso- and hyperendemic areas due to mass drug administration of ivermectin (MDAi). Given the large, contiguous hypo-, meso-, and hyperendemic areas, sustainable elimination of onchocerciasis in sub-Saharan Africa requires delineation of geographic boundaries for parasite transmission zones, so that programs can consider the risk of parasite re-introduction through vector or human migration from areas with ongoing transmission when making decisions to stop MDAi. We propose that transmission zone boundaries can be delineated by characterising the parasite genetic population structure within and between potential zones. We analysed whole mitochondrial genome sequences of 189 O. volvulus adults to determine the pattern of genetic similarity across three West African countries: Ghana, Mali, and Côte d'Ivoire. Population genetic structure indicates that parasites from villages near the Pru, Daka, and Black Volta rivers in central Ghana belong to one parasite population, indicating that the assumption that river basins constitute individual transmission zones is not supported by the data. Parasites from Mali and Côte d'Ivoire are genetically distinct from those from Ghana. This research provides the basis for developing tools for elimination programs to delineate transmission zones, to estimate the risk of parasite re-introduction via vector or human movement when intervention is stopped in one area while transmission is ongoing in others, to identify the origin of infections detected post-treatment cessation, and to investigate whether persisting prevalence despite ongoing interventions in one area is due to parasites imported from others.

Strongyloides questions-a research agenda for the future.

The Strongyloides genus of parasitic nematodes have a fascinating life cycle and biology, but are also important pathogens of people and a World Health Organization-defined neglected tropical disease. Here, a community of Strongyloides researchers have posed thirteen major questions about Strongyloides biology and infection that sets a Strongyloides research agenda for the future. This article is part of the Theo Murphy meeting issue 'Strongyloides: omics to worm-free populations'.

Procedure for Handling and Storage of Onchocerca volvulus Microfilariae Obtained from Skin Snips for Downstream Genetic Work.

WHO and endemic countries target elimination of transmission of Onchocerca volvulus, the parasite causing onchocerciasis. Population genetic analysis of O. volvulus may provide data to improve the evidence base for decisions on when, where, and for how long to deploy which interventions and post-intervention surveillance to achieve elimination. Development of necessary methods and tools requires parasites suitable for genetic analysis. Based on our experience with microfilariae obtained from different collaborators, we developed a microfilariae transfer procedure for large-scale studies in the Democratic Republic of Congo (DRC) comparing safety and efficacy of ivermectin, the mainstay of current onchocerciasis elimination strategies, and moxidectin, a new drug. This procedure is designed to increase the percentage of microfilariae in skin snips suitable for genetic analysis, improve assignment to metadata, and minimize time and materials needed by the researchers collecting the microfilariae. Among 664 microfilariae from South Sudan, 35.7% and 39.5% failed the mitochondrial and nuclear qPCR assay. Among the 576 microfilariae from DRC, 16.0% and 16.7% failed these assays, respectively. This difference may not only be related to the microfilariae transfer procedure but also to other factors, notably the ethanol concentration in the tubes in which microfilariae were stored (64% vs. ≥75%).

Assessing Onchocerca volvulus Intensity of Infection and Genetic Diversity Using Mitochondrial Genome Sequencing of Single Microfilariae Obtained before and after Ivermectin Treatment.

Onchocerciasis is a neglected tropical disease targeted for elimination using ivermectin mass administration. Ivermectin kills the microfilariae and temporarily arrests microfilariae production by the macrofilariae. We genotyped 436 microfilariae from 10 people each in Ituri, Democratic Republic of the Congo (DRC), and Maridi County, South Sudan, collected before and 4-5 months after ivermectin treatment. Population genetic analyses identified 52 and 103 mitochondrial DNA haplotypes among the microfilariae from DRC and South Sudan, respectively, with few haplotypes shared between people. The percentage of genotype-based correct assignment to person within DRC was ~88% and within South Sudan ~64%. Rarefaction and extrapolation analysis showed that the genetic diversity in DRC, and even more so in South Sudan, was captured incompletely. The results indicate that the per-person adult worm burden is likely higher in South Sudan than DRC. Analyses of haplotype data from a subsample (n = 4) did not discriminate genetically between pre- and post-treatment microfilariae, confirming that post-treatment microfilariae are not the result of new infections. With appropriate sampling, mitochondrial haplotype analysis could help monitor changes in the number of macrofilariae in a population as a result of treatment, identify cases of potential treatment failure, and detect new infections as an indicator of continuing transmission.

Geospatial modeling of pre-intervention nodule prevalence of Onchocerca volvulus in Ethiopia as an aid to onchocerciasis elimination.

BACKGROUND: Onchocerciasis is a neglected tropical filarial disease transmitted by the bites of blackflies, causing blindness and severe skin lesions. The change in focus for onchocerciasis management from control to elimination requires thorough mapping of pre-control endemicity to identify areas requiring interventions and to monitor progress. Onchocerca volvulus nodule prevalence in sub-Saharan Africa is spatially continuous and heterogeneous, and highly endemic areas may contribute to transmission in areas of low endemicity or vice-versa. Ethiopia is one such onchocerciasis-endemic country with heterogeneous O. volvulus nodule prevalence, and many districts are still unmapped despite their potential for onchocerciasis transmission. METHODOLOGY/PRINCIPLE FINDINGS: A Bayesian geostatistical model was fitted for retrospective pre-intervention nodule prevalence data collected from 916 unique sites and 35,077 people across Ethiopia. We used multiple environmental, socio-demographic, and climate variables to estimate the pre-intervention prevalence of O. volvulus nodules across Ethiopia and to explore their relationship with prevalence. Prevalence was high in southern and northwestern Ethiopia and low in Ethiopia's central and eastern parts. Distance to the nearest river (RR: 0.9850, 95% BCI: 0.9751-0.995), precipitation seasonality (RR: 0.9837, 95% BCI: 0.9681-0.9995), and flow accumulation (RR: 0.9586, 95% BCI: 0.9321-0.9816) were negatively associated with O. volvulus nodule prevalence, while soil moisture (RR: 1.0218, 95% BCI: 1.0135-1.0302) was positively associated. The model estimated the number of pre-intervention cases of O. volvulus nodules in Ethiopia to be around 6.48 million (95% BCI: 3.53-13.04 million). CONCLUSIONS/SIGNIFICANCE: Nodule prevalence distribution was correlated with habitat suitability for vector breeding and associated biting behavior. The modeled pre-intervention prevalence can be used as a guide for determining priorities for elimination mapping in regions of Ethiopia that are currently unmapped, most of which have comparatively low infection prevalence.

Genetic epidemiology of lymphatic filariasis in American Samoa after mass drug administration.

Over 892 million people in 48 countries are at risk of infection by nematodes that cause lymphatic filariasis. As part of the Global Programme to Eliminate Lymphatic Filariasis, mass drug administration is distributed to communities until surveillance indicates infection rates are below target prevalence thresholds. In some countries, including American Samoa, lymphatic filariasis transmission persists despite years of mass drug administration and/or has resurged after cessation. Nothing is known about the population genetics of Wuchereria bancrofti worms in Polynesia, or whether local transmission is persisting and/or increasing due to inadequate mass drug administration coverage, expansion from residual hotspots, reintroduction from elsewhere, or a combination. We extracted DNA from microfilariae on blood slides collected during prevalence surveys in 2014 and 2016, comprising 31 pools of five microfilariae from 22 persons living in eight villages. We sequenced 1104 bp across three mitochondrial markers (ND4, COI, CYTB). We quantified parasite genetic differentiation using variant calls and estimated haplotypes using principal components analysis, F-statistics, and haplotype networks. Of the variants called, all but eight were shared across the main island of Tutuila, and three of those were from a previously described hotspot village, Fagali'i. Genotypic data did not support population genetic structure among regions or villages in 2016, although differences were observed between worms collected in Fagali'i in 2014 and those from 2016. Because estimated haplotype frequency varied between villages, these statistics suggested genetic differentiation, but were not consistent among villages. Finally, haplotype networks demonstrated American Samoan sequence clusters were related to previously published sequences from Papua New Guinea. These are, to our knowledge, the first reports of W. bancrofti genetic variation in Polynesia. The resurgent parasites circulating on the main island of American Samoa represent a single population. This study is the first step towards investigating how parasite population structure might inform strategies to manage resurgence and elimination of lymphatic filariasis.

Lymphatic filariasis epidemiology in Samoa in 2018: Geographic clustering and higher antigen prevalence in older age groups.

BACKGROUND: Samoa conducted eight nationwide rounds of mass drug administration (MDA) for lymphatic filariasis (LF) between 1999 and 2011, and two targeted rounds in 2015 and 2017 in North West Upolu (NWU), one of three evaluation units (EUs). Transmission Assessment Surveys (TAS) were conducted in 2013 (failed in NWU) and 2017 (all three EUs failed). In 2018, Samoa was the first in the world to distribute nationwide triple-drug MDA using ivermectin, diethylcarbamazine, and albendazole. Surveillance and Monitoring to Eliminate LF and Scabies from Samoa (SaMELFS Samoa) is an operational research program designed to evaluate the effectiveness of triple-drug MDA on LF transmission and scabies prevalence in Samoa, and to compare the usefulness of different indicators of LF transmission. This paper reports results from the 2018 baseline survey and aims to i) investigate antigen (Ag) prevalence and spatial epidemiology, including geographic clustering; ii) compare Ag prevalence between two different age groups (5-9 years versus ≥10 years) as indicators of areas of ongoing transmission; and iii) assess the prevalence of limb lymphedema in those aged ≥15 years. METHODS: A community-based cluster survey was conducted in 30 randomly selected and five purposively selected clusters (primary sampling units, PSUs), each comprising one or two villages. Participants were recruited through household surveys (age ≥5 years) and convenience surveys (age 5-9 years). Alere Filariasis Test Strips (FTS) were used to detect Ag, and prevalence was adjusted for survey design and standardized for age and gender. Adjusted Ag prevalence was estimated for each age group (5-9, ≥10, and all ages ≥5 years) for random and purposive PSUs, and by region. Intraclass correlation (ICC) was used to quantify clustering at regions, PSUs, and households. RESULTS: A total of 3940 persons were included (1942 children aged 5-9 years, 1998 persons aged ≥10 years). Adjusted Ag prevalence in all ages ≥5 years in randomly and purposively selected PSUs were 4.0% (95% CI 2.8-5.6%) and 10.0% (95% CI 7.4-13.4%), respectively. In random PSUs, Ag prevalence was lower in those aged 5-9 years (1.3%, 95% CI 0.8-2.1%) than ≥10 years (4.7%, 95% CI 3.1-7.0%), and poorly correlated at the PSU level (R-square = 0.1459). Adjusted Ag prevalence in PSUs ranged from 0% to 10.3% (95% CI 5.9-17.6%) in randomly selected and 3.8% (95% CI 1.3-10.8%) to 20.0% (95% CI 15.3-25.8%) in purposively selected PSUs. ICC for Ag-positive individuals was higher at households (0.46) compared to PSUs (0.18) and regions (0.01). CONCLUSIONS: Our study confirmed ongoing transmission of LF in Samoa, in accordance with the 2017 TAS results. Ag prevalence varied significantly between PSUs, and there was poor correlation between prevalence in 5-9 year-olds and older ages, who had threefold higher prevalence. Sampling older age groups would provide more accurate estimates of overall prevalence, and be more sensitive for identifying residual hotspots. Higher prevalence in purposively selected PSUs shows local knowledge can help identify at least some hotspots.

Cryptic species diversity in ticks that transmit disease in Australia.

Ticks are important vectors of a broad range of pathogens in Australia. Many tick species are morphologically similar and are therefore difficult to identify using morphology alone, particularly when collected in the larval and nymphal life stages. We report here the application of molecular methods to examine the species diversity of ixodid ticks at two sites in southern New South Wales, Australia. Our taxon sampling included six morphologically characterised adult stage voucher specimens of Ixodes trichosuri, Ixodes tasmani, Ixodes fecialis and Ixodes holocyclus (the paralysis tick) and ~250 field collected specimens that were in the larva or nymph stage and thus not morphologically identifiable. One nuclear and two mitochondrial amplicons were sequenced using a combination of Sanger and Illumina MiSeq sequencing. Phylogenetic relationships were estimated using both maximum likelihood and Bayesian methods. Two clades with strong bootstrap and Bayesian support were observed across trees estimated from each of three markers and from an analysis of the concatenated sequences. One voucher specimen of I. trichosuri was located in one of these clades, while the other I. trichosuri voucher specimen was in a second clade with the remaining three identified species, suggesting these morphologically similar ticks may represent different cryptic species. Unidentified specimens were found across both clades, and molecular divergence of many of these is equal to or greater than that observed between identified species, suggesting additional unidentified species may exist. Further studies are required to understand the taxonomic status of ticks in Australia, and how this species diversity impacts disease risk for livestock, domestic animals, wildlife and humans.

Genomic Epidemiology in Filarial Nematodes: Transforming the Basis for Elimination Program Decisions.

Onchocerciasis and lymphatic filariasis are targeted for elimination, primarily using mass drug administration at the country and community levels. Elimination of transmission is the onchocerciasis target and global elimination as a public health problem is the end point for lymphatic filariasis. Where program duration, treatment coverage, and compliance are sufficiently high, elimination is achievable for both parasites within defined geographic areas. However, transmission has re-emerged after apparent elimination in some areas, and in others has continued despite years of mass drug treatment. A critical question is whether this re-emergence and/or persistence of transmission is due to persistence of local parasites-i.e., the result of insufficient duration or drug coverage, poor parasite response to the drugs, or inadequate methods of assessment and/or criteria for determining when to stop treatment-or due to re-introduction of parasites via human or vector movement from another endemic area. We review recent genetics-based research exploring these questions in Onchocerca volvulus, the filarial nematode that causes onchocerciasis, and Wuchereria bancrofti, the major pathogen for lymphatic filariasis. We focus in particular on the combination of genomic epidemiology and genome-wide associations to delineate transmission zones and distinguish between local and introduced parasites as the source of resurgence or continuing transmission, and to identify genetic markers associated with parasite response to chemotherapy. Our ultimate goal is to assist elimination efforts by developing easy-to-use tools that incorporate genetic information about transmission and drug response for more effective mass drug distribution, surveillance strategies, and decisions on when to stop interventions to improve sustainability of elimination.

Identifying residual transmission of lymphatic filariasis after mass drug administration: Comparing school-based versus community-based surveillance - American Samoa, 2016.

INTRODUCTION: Under the Global Programme to Eliminate Lymphatic Filariasis (LF), American Samoa conducted seven rounds of mass drug administration (MDA) from 2000-2006. The World Health Organization recommends systematic post-MDA surveillance using Transmission Assessment Surveys (TAS) for epidemiological assessment of recent LF transmission. We compared the effectiveness of two survey designs for post-MDA surveillance: a school-based survey of children aged 6-7 years, and a community-based survey targeting people aged ≥8 years. METHODS: In 2016, we conducted a systematic school-based TAS in all elementary schools (N = 29) and a cluster survey in 28 villages on the two main islands of American Samoa. We collected information on demographics and risk factors for infection using electronic questionnaires, and recorded geo-locations of schools and households. Blood samples were collected to test for circulating filarial antigen (CFA) using the Alere Filariasis Test Strip. For those who tested positive, we prepared slides for microscopic examination of microfilaria and provided treatment. Descriptive statistics were performed for questionnaire variables. Data were weighted and adjusted to account for sampling design and sex for both surveys, and for age in the community survey. RESULTS: The school-based TAS (n = 1143) identified nine antigen-positive children and found an overall adjusted CFA prevalence of 0.7% (95% CI: 0.3-1.8). Of the nine positive children, we identified one microfilariaemic 7-year-old child. The community-based survey (n = 2507, 711 households) identified 102 antigen-positive people, and estimated an overall adjusted CFA prevalence of 6.2% (95% CI: 4.5-8.6). Adjusted village-level prevalence ranged from 0-47.1%. CFA prevalence increased with age and was higher in males. Of 86 antigen-positive community members from whom slides were prepared, 22 (25.6%) were microfilaraemic. School-based TAS had limited sensitivity (range 0-23.8%) and negative predictive value (range 25-83.3%) but had high specificity (range 83.3-100%) and positive predictive value (range 0-100%) for identifying villages with ongoing transmission. CONCLUSIONS: American Samoa failed the school-based TAS in 2016, and the community-based survey identified higher than expected numbers of antigen-positive people. School-based TAS was logistically simpler and enabled sampling of a larger proportion of the target population, but the results did not provide a good indication of the overall CFA prevalence in older age groups and was not sensitive at identifying foci of ongoing transmission. The community-based survey, although operationally more challenging, identified antigen-positive individuals of all ages, and foci of high antigen prevalence. Both surveys confirmed recrudescence of LF transmission.

Introduction of Non-Native Pollinators Can Lead to Trans-Continental Movement of Bee-Associated Fungi.

Bees are essential pollinators for many flowering plants, including agriculturally important crops such as apple. As geographic ranges of bees or their host plants change as a result of human activities, we need to identify pathogens that could be transmitted among newly sympatric species to evaluate and anticipate their effects on bee communities. We used PCR screening and DNA sequencing to evaluate exposure to potentially disease-causing microorganisms in a pollinator of apple, the horned mason bee (Osmia cornifrons). We did not detect microsporidia, Wolbachia, or trypanosomes, which are common pathogens of bees, in any of the hundreds of mason bees screened. We did detect both pathogenic and apathogenic (saprophytic) fungal species in the genus Ascosphaera (chalkbrood), an unidentified species of Aspergillus fungus, and a strain of bacteria in the genus Paenibacillus that is probably apathogenic. We detected pathogenic fungal strains in asymptomatic adult bees that therefore may be carriers of disease. We demonstrate that fungi from the genus Ascosphaera have been transported to North America along with the bee from its native range in Japan, and that O. cornifrons is exposed to fungi previously only identified from nests of other related bee species. Further study will be required to quantify pathogenicity and health effects of these different microbial species on O. cornifrons and on closely-related native North American mason bees that may now be exposed to novel pathogens. A global perspective is required for pathogen research as geographic ranges of insects and microorganisms shift due to intentional or accidental introductions.

The bee tree of life: a supermatrix approach to apoid phylogeny and biogeography.

BACKGROUND: Bees are the primary pollinators of angiosperms throughout the world. There are more than 16,000 described species, with broad variation in life history traits such as nesting habitat, diet, and social behavior. Despite their importance as pollinators, the evolution of bee biodiversity is understudied: relationships among the seven families of bees remain controversial, and no empirical global-level reconstruction of historical biogeography has been attempted. Morphological studies have generally suggested that the phylogeny of bees is rooted near the family Colletidae, whereas many molecular studies have suggested a root node near (or within) Melittidae. Previous molecular studies have focused on a relatively small sample of taxa (~150 species) and genes (seven at most). Public databases contain an enormous amount of DNA sequence data that has not been comprehensively analysed in the context of bee evolution. RESULTS: We downloaded, aligned, concatenated, and analysed all available protein-coding nuclear gene DNA sequence data in GenBank as of October, 2011. Our matrix consists of 20 genes, with over 17,000 aligned nucleotide sites, for over 1,300 bee and apoid wasp species, representing over two-thirds of bee genera. Whereas the matrix is large in terms of number of genes and taxa, there is a significant amount of missing data: only ~15% of the matrix is populated with data. The placement of the root as well as relationships between Andrenidae and other bee families remain ambiguous, as several alternative maximum-likelihood estimates fall within the statistically credible set. However, we recover strong bootstrap support for relationships among many families and for their monophyly. Ancestral geographic range reconstruction suggests a likely origin of bees in the southern hemisphere, with Melittidae ancestrally located within Africa, and Halictidae, Colletidae, and Apidae within the New World. CONCLUSIONS: Our study affirms the monophyly of each bee family, sister-taxa relationships between Apidae and Megachilidae (the 'long-tongued bees'), between Colletidae and Stenotritidae, and between Colletidae + Stenotritidae and Halictidae. Our analyses reject a Colletidae-basal hypothesis for family-level relationships and instead support Melittidae as sister to the remaining bees. Southern hemisphere vicariance likely played an important role in early diversification within many bee families.

Targeted enrichment: maximizing orthologous gene comparisons across deep evolutionary time.

Estimated phylogenies of evolutionarily diverse taxa will be well supported and more likely to be historically accurate when the analysis contains large amounts of data-many genes sequenced across many taxa. Inferring such phylogenies for non-model organisms is challenging given limited resources for whole-genome sequencing. We take advantage of genomic data from a single species to test the limits of hybridization-based enrichment of hundreds of exons across frog species that diverged up to 250 million years ago. Enrichment success for a given species depends greatly on the divergence time between it and the reference species, and the resulting alignment contains a significant proportion of missing data. However, our alignment generates a well-supported phylogeny of frogs, suggesting that this technique is a practical solution towards resolving relationships across deep evolutionary time.

Historical changes in northeastern US bee pollinators related to shared ecological traits.

Pollinators such as bees are essential to the functioning of terrestrial ecosystems. However, despite concerns about a global pollinator crisis, long-term data on the status of bee species are limited. We present a long-term study of relative rates of change for an entire regional bee fauna in the northeastern United States, based on >30,000 museum records representing 438 species. Over a 140-y period, aggregate native species richness weakly decreased, but richness declines were significant only for the genus Bombus. Of 187 native species analyzed individually, only three declined steeply, all of these in the genus Bombus. However, there were large shifts in community composition, as indicated by 56% of species showing significant changes in relative abundance over time. Traits associated with a declining relative abundance include small dietary and phenological breadth and large body size. In addition, species with lower latitudinal range boundaries are increasing in relative abundance, a finding that may represent a response to climate change. We show that despite marked increases in human population density and large changes in anthropogenic land use, aggregate native species richness declines were modest outside of the genus Bombus. At the same time, we find that certain ecological traits are associated with declines in relative abundance. These results should help target conservation efforts focused on maintaining native bee abundance and diversity and therefore the important ecosystems services that they provide.

Multiple origins of kdr-type resistance in the house fly, Musca domestica.

Insecticide resistance is a model phenotype that can be used to investigate evolutionary processes underlying the spread of alleles across a global landscape, while offering valuable insights into solving the problems that resistant pests present to human health and agriculture. Pyrethroids are one of the most widely used classes of insecticides world-wide and they exert their toxic effects through interactions with the voltage-sensitive sodium channel (Vssc). Specific mutations in Vssc (kdr, kdr-his and super-kdr) are known to cause resistance to pyrethroid insecticides in house flies. In order to determine the number of evolutionary origins of kdr, kdr-his and super-kdr, we sequenced a region of Vssc from house flies collected in the USA, Turkey and China. Our phylogenetic analysis of Vssc unequivocally supports the hypothesis of multiple independent origins of kdr, super-kdr and kdr-his on an unprecedented geographic scale. The implications of these evolutionary processes on pest management are discussed.

Rare gene capture in predominantly androgenetic species.

The long-term persistence of completely asexual species is unexpected. Although asexuality has short-term evolutionary advantages, a lack of genetic recombination leads to the accumulation over time of deleterious mutations. The loss of individual fitness as a result of accumulated deleterious mutations is expected to lead to reduced population fitness and possible lineage extinction. Persistent lineages of asexual, all-female clones (parthenogenetic and gynogenetic species) avoid the negative effects of asexual reproduction through the production of rare males, or otherwise exhibit some degree of genetic recombination. Another form of asexuality, known as androgenesis, results in offspring that are clones of the male parent. Several species of the Asian clam genus Corbicula reproduce via androgenesis. We compared gene trees of mitochondrial and nuclear loci from multiple sexual and androgenetic species across the global distribution of Corbicula to test the hypothesis of long-term clonality of the androgenetic species. Our results indicate that low levels of genetic capture of maternal nuclear DNA from other species occur within otherwise androgenetic lineages of Corbicula. The rare capture of genetic material from other species may allow androgenetic lineages of Corbicula to mitigate the effects of deleterious mutation accumulation and increase potentially adaptive variation. Models comparing the relative advantages and disadvantages of sexual and asexual reproduction should consider the possibility of rare genetic recombination, because such events seem to be nearly ubiquitous among otherwise asexual species.

When trees grow too long: investigating the causes of highly inaccurate bayesian branch-length estimates.

A surprising number of recent Bayesian phylogenetic analyses contain branch-length estimates that are several orders of magnitude longer than corresponding maximum-likelihood estimates. The levels of divergence implied by such branch lengths are unreasonable for studies using biological data and are known to be false for studies using simulated data. We conducted additional Bayesian analyses and studied approximate-posterior surfaces to investigate the causes underlying these large errors. We manipulated the starting parameter values of the Markov chain Monte Carlo (MCMC) analyses, the moves used by the MCMC analyses, and the prior-probability distribution on branch lengths. We demonstrate that inaccurate branch-length estimates result from either 1) poor mixing of MCMC chains or 2) posterior distributions with excessive weight at long tree lengths. Both effects are caused by a rapid increase in the volume of branch-length space as branches become longer. In the former case, both an MCMC move that scales all branch lengths in the tree simultaneously and the use of overdispersed starting branch lengths allow the chain to accurately sample the posterior distribution and should be used in Bayesian analyses of phylogeny. In the latter case, branch-length priors can have strong effects on resulting inferences and should be carefully chosen to reflect biological expectations. We provide a formula to calculate an exponential rate parameter for the branch-length prior that should eliminate inference of biased branch lengths in many cases. In any phylogenetic analysis, the biological plausibility of branch-length output must be carefully considered.